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ABSTRACT 

We discuss the acceleration and escape of secondary particles, especially 
positrons produced by hadronic interactions in a supernova remnant (SNR) 
shock. During the shock acceleration, protons would interact with ambient gas 
and produce charged secondary particles, which would also be accelerated in 
a SNR and injected into the interstellar medium as cosmic-rays (CRs). Some 
previous studies proposed that the positron excess observed by PAMELA can 
be explained with this process. We calculate the energy spectra of CR protons 
and secondary CR positrons escaping a SNR into the interstellar medium in a 
time-dependent manner. We show that, on the contrary to the results presented 
previously, the observed spectra of secondary CR particles generated in SNRs 
would be softer than those of primary CR particles, which means that the rise 
in the positron fraction cannot be reproduced by this model. This is because 
protons with higher energy would be trapped around the shock only for a shorter 
time compared to those with lower energy when the maximum CR energy is 
escape-limited, which results in a smaller number of secondary CR particles with 
higher energy. 

Subject headings: acceleration of particles - cosmic rays - supernova remnants - 
shock waves 



1. Introduction 



Recently the PAMELA collaboration reported that the cosmic-ray (CR) positron frac- 
tion e + / (e + + e~) is rising with energy in the range of ~ 7 — lOOGeV (Adriani et al. 2009), 
which cannot be reproduced when considering only the positron production during the prop- 
agation of CR protons in the Galaxy. After that, Fermi Large Area Telescope (LAT) 
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observed the electron spectrum and positron spectrum separately making use of Earth's 
magnetic field, and confirmed this trend in the energy range of 20-200GeV (Ackermann et 
al. 2012). In addition, the cosmic-ray electrons plus positrons flux has been measured by 
Fermi/H.E.S.S./ATIC/PPB-BETS and they are also shown to have the excess in compari- 
son with the standard CR propagation model (Abdo et al. 2009a; Ackermann et al. 2010; 
Aharonian et al. 2008,2009; Chang et al. 2008; Torii et al. 2008). These observations likely 
suggest a new kind of source of CR electrons and positrons, such as nearby pulsars (Profumo 
2008; Malyshev et al. 2009; Grasso et al. 2009; Hooper et al. 2009; Yuksel et al. 2009; 
Kawanaka et al. 2010; Heyl et al. 2010; Kisaka & Kawanaka 2012), microquasars (Heinz 
& Sunyaev 2002), gamma-ray bursts (Ioka 2010; Calvez & Kusenko 2010), and dark matter 
annihilations/decays (see Bergstrom et al. 2008; Arkani-Hamed et al. 2009; Maede et al. 
2010; Nardi et al. 2009; Grasso et al. 2009 and references therein). 

Supernova remnants (SNRs) are also promising candidates of CR positron sources (Fu- 
jita et al. 2009; Shaviv et al. 2009; Hu et al. 2009; Blasi 2009; Biermann et al. 2009; 
Lee et al. 2011). It is widely believed that Galactic CR nuclei and electrons are acceler- 
ated in SNRs by the diffusive shock acceleration (DSA) mechanism (Bell 1978; Blandford & 
Ostriker 1978). The theory of DSA (for reviews, see Blandford & Eichler 1987; Malkov & 
Drury 2001) can naturally derive the power-law spectrum of particles accelerated in the SNR 
shock. The TeV gamma-ray detections from the shell of young SNRs by H.E.S.S. (Aharonian 
et al. 2004, 2005) and the GeV gamma-ray detections from middle-aged SNRs interacting 
with molecular clouds by Fermi & AGILE (Abdo et al. 2009b, 2010a, 2010b, 2010c; Tavani 
et al. 2010; Giuliani et al. 2010) can be interpreted as the results of hadronic interactions 
between CR protons accelerated in SNRs and ambient protons. 

Blasi (2009) has proposed the model of positron production which considers the hadronic 
interactions during the proton acceleration in the SNR shock. In this model the acceleration 
of primary CR protons and the production and acceleration of secondary positrons take place 
in the same region. According to this study the resulting positron spectrum would be harder 
than the primary electron spectrum (locally having the same index as the primary protons) 
and therefore the positron fraction would rise with energy as the data reported by PAMELA. 
Because this model can also predict the production of antiprotons and secondary nuclei (B, 
Ti etc.), it is suggested that antiproton-to-proton ratio (p/p), boron-to-carbon ratio (B/C) 
and titanium-to-iron ratio (Ti/Fe) would show the excess above ~ lOOGeV per nucleon 
(Blasi & Serpico 2009; Mertsch & Sarkar 2009; Ahlers et al. 2009; see also Kachelriess et 
al. 2011, who argued that this process cannot reproduce the rising of the positron fraction 
using Monte Carlo simulations). These predictions in CR spectral features would be tested 
by the ongoing observations by AMS-02, together with the more precise data of the positron 
fraction which will be also provided after its measurement. 
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In order to predict the spectrum of these CR components observed at the Earth, we 
should calculate their outgoing flux from SNRs which generally should vary with time accord- 
ing to the dynamical evolution of SNRs. The previous studies shown above, however, only 
estimate the spectrum of CRs advected to the downstream of the SNR shock with totally 
time-independent picture. In this study, we discuss the secondary CR positron production 
process inside the SNR taking into account the time-dependent escape of them as well as 
that of primary CR protons. The escape of CR particles from SNRs and the observational 
feature expected from that have been investigated by several authors (Ptuskin & Zirakashvili 
2005; Reville et al. 2009; Ohira et al. 2010, 2011a, 2011b; Caprioli et al. 2010; Drury 2011; 
Kawanaka et al. 2011). According to these studies, the particles accelerated in the shock 
would escape the shock in a time-dependent manner - the higher energy particles can escape 
the SNR earlier. As the shock slows down and the magnetic field decays, the lower energy 
particles can leave the shock gradually. This process can make the CR spectrum in the in- 
terstellar medium softer than that inside the SNRs (Ohira et al. 2010; Caprioli et al. 2010), 
which can give the explanation for the observed CR spectral index being softer than that 
predicted from DSA theory (~ 2) and the spectral breaks found in the GeV emissions from 
the middle-aged SNRs interacting with molecular clouds (Ohira et al. 2011). Within this 
picture, secondary particles produced and accelerated in the SNR shock would also escape 
the source in a time-dependent manner. Therefore, in order to predict the observed spectra 
of those secondary particles, we should calculate their time-dependent escape from SNRs, as 
well as their production and acceleration. 

This paper is structured as follows: In Section 2, we describe the model of time- 
dependent escape of primary and secondary cosmic-ray particles from SNRs and show the 
transport equations of them. In Section 3, we solve the transport equation for secondary 
cosmic-ray positrons, and predict their energy spectrum. We discuss our results in Section 
4, and Section 5 is devoted to the summary or this work. 

2. Model 

According to the standard theory of DSA, the scattering of particles is due to the 
turbulent magnetic field around the shock. Especially, in the upstream region, the turbulence 
is considered to be excited by the accelerated particles themselves (Bell 1978). Then, in the 
far upstream region which is sufficiently distant from the shock front, because the particle flux 
should be small, the turbulence would be so weak that the particles cannot be backscattered 
to the shock and escape the SNR. In other words, in the CR escape scenario, we should 
impose the escape boundary in the upstream with a finite distance (hereafter, /) from the 
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shock front. Letting D{p) be the diffusion coefficient for a particle with momentum of p and 
u sh be the velocity of the shock front, the escape condition for a particle can be described as 



where the left hand side is the characteristic diffusion length of a particle with momentum p. 
The escape boundary x = —I (here x < is the upstream region) is the point beyond which 
particle-scattering Alfven waves are absent and therefore the mean free path of particles is 
much larger than the size of a SNR. Zirakashvili & Ptuskin (2008) argued that the distance 
between the escape boundary and the shock front should be given roughly by the radius of 
the remnant ~ R s h- In fact, because of the spherical geometry of a SNR shock, the particles 
which have diffused out from the shock to the distance comparable with the shock radius 
would not be able to cross the shock again and therefore escape the SNR into the interstellar 
medium (ISM). 



2.1. Escape of Primary Cosmic- Ray Protons 



First let us discuss the escape flux of cosmic-ray protons as well as the maximum at- 
tainable energy of them as functions of time following the procedure presented in Caprioli 
et al. (2009) and Ohira et al. (2010). The transport equation for primary protons in the 
stationary case in the shock rest frame is 

pdu df 
3 dx dp 

where fi(x,p) is the distribution function of protons and Qi is the source term. Here u(x) 
is the fluid velocity. The fluid velocity u(x) is given by 

■ui (x < 0) 



u{x) dfi = d_ 
dx dx 



1 + Qi, 



(2) 



u(x) 



(3) 



u 2 (x > 0), 

where U\ and u 2 are constants. Especially U\ us equal to the speed of the shock front in 
the laboratory frame. When considering the time- dependent escape scenario, we use the 
boundary condition f\(x = —l,p) = 0, and then the solution of this transport equation is 
given by 



fi(x,p) = 



1 — exp(— uil/D) 



cxp 



U\X 



cxp 



D(p) 



where /i,o(p) is determined from the junction condition at x — 0, 



/i,o(p) = K(t)exp 



"7 



dlnp' 



Pmin 1 — eXp 



Ull 

' D(p') 



(4) 



(5) 
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where 7 = 3ui/(ui — U2) and K(t) is the time-dependent normalization factor which will be 



determined belowo Here we can derive the escape flux at the boundary x = — I as 



i(p) 



Mi /1 



* ~ D{P) lx~ 



«i/i,o(p) 



exp(uil / D) 



(6) 



and this value is negative which means that the flux is directed to the upstream (i.e. outward 
of the shock). We assume the Bohm-type diffusion in which the mean free path of a particle 
is proportional to its Larmor radius. Then the diffusion coefficient for particles with a 
momentum p and a charge e can be expressed as 



D(p) = % 



2 

c p 
3eB' 



(7) 



where B is the magnetic field strength and rj g is the gyrofactor which should be equal to 
unity for the Bohm limit. In this case the absolute value of this escape flux —<f>x(p) would 
attain its maximum value when the momentum is 



P 



jDc 



Pi,: 



where Do = D(p)/p, and then the maximum value of the escape flux is 

- MPl,m) = J , (9) 

(Caprioli et al. 2009). As /i,o(p) oc p -7 for p p x ,m an d /i,o oc exp(— p/pi tm ) for p l m ^ p, 
this quantity p\^ m plays the role of maximum momentum of the primary accelerated particles, 
i.e. the accelerated protons with this momentum would escape the SNR into the ISM most 
efficiently According to Ohira et al. (2010), the distribution function of this escape flux can 
be approximated as 



- 01 (?) 



Ml/l,o(Pl,r 
Ulfl,o(pi,r 

e~f - 1 



cxp 



hip — lnpi t 



a 



yGvcr ■ 5 (hip - lnpi, m ), 



(10) 



^^In the realistic situation, the minimum momentum of the integral appearing in Eq.([5]) p mnl should be 
determined by taking into account the particle injection process and should vary with time in general. 
Because the variation of p m i n can be regarded as the change of the overall normalization factor of the 
distribution function /i,o(p), in the following discussion we assume that K(t) varies with time instead of 
Pmin which is fixed as a constant. This assumption does not affect the results presented below. 
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where a 2 = 2(1 — e 7 )/7- 

Because D which appears in the expression of pi j7n is proportional to rjgB -1 , the max- 
imum momentum of escaping particles pi >m depends on the evolution of the amplification 
and decay of the magnetic field and turbulence around the shock (Ptuskin & Zirakashivili 
2003, 2005; Yan et al. 2012). Here we assume that, during the Sedov phase, the maximum 
momentum of escaping primary particles decreases as oc t~ a where a is a parameter which 
characterizes the evolution of pi jTO . Note that, as is obvious from Eq. (|S]), this assumption is 
equivalent with the assumption that the combination of Uil/D is proportional to t~ a . Now 
we consider the Sedov phase and assume that I is proportional to the shock radius i? s h, the 
coefficient D$ evolves with time being proportional to t Q-1//5 . This time dependence can be 
interpreted also in the following way; the maximum momentum of escaping particles should 
be determined by the condition that the timescale for a particle to escape the source with 
size of ~ i? s h, t csc ~ -Rg^Z-DQo), is comparable with the acceleration timescale in the SNR 
shock, t acc ~ D{jS)/u s \y. Considering the time dependence of -Rsh and w s h in the Sedov phase, 
we have the time dependence of the maximum momentum which can escape the source as 
p m3X oc Dq 1 ^ 1 / 5 . Therefore, when imposing that p max oc t~ a , D should be proportional to 

On the other hand, the normalization factor of the escape flux would evolve with time 
approximately as oc uiK(t)pi'J ri . Considering that the energy spectrum of escaping particles 
is described as 

dNesc,l _ , d2 A 2 



dedt 



oc47ri4-4V0iCp), (11) 



where e = ■>/ (cp) 2 + (m p c 2 ) 2 is the energy of a particle, and that R s h and w s h evolve as 
oc t 2//5 and oc £~ 3 / 5 in the Sedov phase, respectively, we can derive the power-law index of 
the time-integrated energy spectrum of primary cosmic-ray protons s esCj i as 

SesCjl ~ ( 7 - 2) + 1 (h + , (12) 

where we use Eq. ( TlOl) and assume that K(t) is proportional to t 13 . If we redefine as s = 7 — 2 
and /3 = /3 + 6/5, this formula is identical to Eq. (28) in Ohira et al. (2010). In the standard 
theory of non-relativistic DSA, the index 7 should be equal to 4 and therefore the energy 
spectral index of accelerated particles should be 2. However, if the normalization factor 
of CR particle flux increases with time (i.e. if /5 > 0), the spectral index of escaping CR 
particles would be softer than 2, which is just what has been implied from the observations 
(Ohira et al. 2010). 
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2.2. Transport Equation of Secondary Cosmic-Ray Particles 

During the acceleration of protons in the SNR shocks, secondary particles such as 
positrons, antiprotons are produced through hadronic interactions. In order to calculate 
their spectrum, we should solve the transport equation with a source term which is related 
to the primary proton distribution function derived above: 



u\x) 



,dh 
dx 



d_ 

dx 



. pdudf 2 ( , 
6 ax op 



(13) 



Here f 2 (x,p) is the distribution function of secondary particles and Q 2 (x,p) is the source 
term which can be determined from Eq.fll]) as 



Q2(x,p) 



Q2M 



cxp 



-exp 



[ mi 

1 CXP L D(p p )\ 

(x > 0), 



(14) 



where p p is the momentum of a secondary particle generated from a primary particle with 
a momentum p. Generally p and p p can be approximately related linearly: p w £jp p . In 



the case of positron production, £ e 



0.05, and in the case of p, £~ ~ 0.17. Here we 



assume D(p) oc p (Bohm diffusion), so we can replace D(p p ) by D(p)/£i. The momentum 
dependences Q2, u {p) and Qz^ip) can be expressed as 



Q2,u{p) 
Q2AP) 



cn. 



gas. 11 



rQ2,u(p), 



1 j 

dp 



(15) 
(16) 



where N\(p) = 4iTp 2 fi(0,p) is the number density of primary CR particles per unit momen- 
tum, dtjpp/ dp is the differential cross section of i particle generations via pp interactions, and 
Pth is the threshold momentum of primary protons for % particle generations. 

The solution of Eq.( lT3"|) should satisfy the following boundary conditions: 



11 



m 



iv 



lim f 2 ,u 

x—t—O 



lim fo , 



D{p) 



dx 



-^-ujp — 



x=0 



lim f 2u = 0, 

x— >— I 



lim f 2 

x— >oo 



< 00. 



Here the condition (ii) comes from the integration of Eq. ffTB"]) across the shock front, and it 
yields the differential equation with respect of p. 
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3. Energy Spectrum of Escaping Secondary Positrons 



Under the boundary conditions shown above, we can solve the transport equation ana- 
lytically as 



f2, u = Vi(p)exp — + Y 2 (p)exp 



D 



+ Y 3 (p)x + Y±(p), 



h,d = f2{X = 0,P) + X, 



u 2 



where 



Yi(p) 
Y 2 (p) 
Y 3 (p) 
Y^P) 



1 — exp(— u\l/D) 
D Q 2 ,M 



' , , n , D Q 2 , u (p) I Q2M 

J2(}J,P) ^—TZ TT + 



«i &(1 - &) u i exp(^«iZ/D) - 1 



1 



i <92,u(p) 



■ui exp(£iUil / D) — 1 ' 
1 



/ 2 (0,p) ( 1 
4 



1 — exp(— uil/D) 
I (D Q2M I Q2M 



1 - exp(-«iZ/D) V M i 6 - C 2 u i exp(^«iZ/D) - 1 
£ Q 2)tt (p) 1 



/ 2 (0,p) 
G(p) 



= -7 



u? 1 - exp(-&itiZ/£>) & - Q ' 

" P dp' if (p) Q2,U(P'W) 



H(p') 



1 



1 — exp(— uil/D) 



Ul 

D 1 



+ 



Z 



#(p) 



+— 

Mi 

exp 



1 + 



-7 



6 



Mi & - C 2 exp(^uiZ/D) - 1 
r 2 D 



1 - & 1 - exp(-^uiZ/D) 
p dp' 1 



Ul 



'Pmin P' l-exp(-«iZ/D). 
and so the escape flux of the secondary particles at x — — I is 

9f2,u 



4>2 (p) = «l/2,uU=-I - #(p) 



x=—l 



Ul 



exp(u\l/D) — 1 



f 2 (o,p) - —- — - + 



<?2,u 



u\ ii - Q mi exp(^uiZ/D) - 1 



Q2,v 



M^l-^exp^MiZ/L*)-!' 



(17) 
(18) 

(19) 
(20) 
(21) 



(22) 
(23) 



(24) 
(25) 



(26) 
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and in order to derive the spectrum of the CR particles injected from a SNR into the 
interstellar medium, we should integrate this escaping CR flux with time. 

In the following calculations, we assume the explosion energy of a supernova as = 
10 51 erg, the number density of the interstellar matter as n = 1cm -3 , and the ejecta mass as 
M e j = IMq. With this environmental parameters, the SNR shock radius at the beginning 
of the Sedov phase, Rs, and the SNR age at that time ts can be expressed as 

R ° = 213pc {it) GsfO • (27) 
is = 209yr (lo^J [m-J (w) ' (28) 

In addition, we perform the calculations assuming a = 2.6 and 

cp lim = 10 6 - 5 GeV x l-j , (29) 

which is also adopted in Gabici et al. (2009) and Ohira et al. (2010) according to the 
hypothesis that SNRs are responsible for the observed CRs with energies from ~ lGeV up 
to the knee energy (~ 10 15 ' 5 GeV). In this case, the lowest energy CR particles (e ~ lGeV) 
would escape the SNR at the end of the Sedov phase t = 10 2,5 ts, when the shock radius R^ 
becomes 10-Rs. Note that the magnetic field strength at the beginning of the Sedov phase 
should be amplified up to 

Bs — 174^G (3^) (rf) (^) (^)^(^)^ (30) 



In calculating the positron production from hadronic interactions in the SNR by Eq. (TT5]) . 
we used the parametrization for the cross section provided by Kamae et al. (2006). Also, 
we assume that r = 4 and 7 = 4, which is true in the case of non-relativistic strong shock, 
which is realized in a typical SNR. 



3.1. Escape Flux 

In Figure 1, we show the escape fluxes of primary CR protons and secondary CR 
positrons as functions of momentum. In the latter one, we can see two peaks at differ- 
ent momentums. These peaks come from two dominant terms in Eq. (T26l) : 

, / A _ "1/2(0, p) 

-<fa,A[P) = ( ;/n x 7> 31 

exp(wi</ D) — 1 



- £ e +) exp(^ e +ni//D) - 1 ; 
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where 2 ,aG°) contributes to the peak at a higher momentum and <p2,B{p) to the peak at 
a lower momentum. Judging from their functional forms and the figure, <p2,A has a similar 
form to the escape flux of primary protons 0i, and so this component can be regarded as 
the positrons accelerated at the shock in the same manner as protons. On the other hand, 
<p2,B has a peak with momentum of ~ 0.05pi )inax ~ £ e +Pi,max- Therefore, we can say that the 
positrons in this component has been produced mainly by protons with momentum around 
Pi.max- Note that, though in Fig. 1 the peak of 2 ,s seems to be higher than that of 2 , b, the 
escape flux per unit energy bin is proportional to p 2 x 2 (p), and then the main contribution 
to the resulting energy spectrum of escaping positrons comes from (p2,A, not from 2 ,s (see 
below) . 

In Figure 2 we show the energy distribution functions of the escape flux (oc p 2 <p{p)) of 
primary protons (upper panel) and secondary positrons (lower panel) at t = I0t s , 30t s , 100tg, 
and 300ts- In these calculations, we neglect the energy loss of positrons during acceleration 
due to synchrotron radiation. If we take into account this effect, positrons can escape the 
SNR only when their escape timescale is shorter than their cooling timescale: 

R\ 9m 4 c 7 
tesc tc °° l ~ 4ABfe? (33) 

where m e and e e + are the mass and the energy of a positron, respectively, and Bd = AB 
is the magnetic field strength in the downstream region. Ohira et al. (2011) showed that 
the highest energy of electrons which can escape the SNR is E m e ~ 50TeV. Therefore, we 
cannot observe the positrons which are more energetic than this E m}B . As long as discussing 
positrons with energy lower than this value, the results shown here describe the phenomena 
well. 

We can see that, (1) the energy distribution function of the secondary CR positron flux 
has a peak at a certain energy in each time step, and that (2) the peak energy decreases 
with time in the same rate as that of primary protons. The fact (1) means that, as in the 
case of primary CR protons, there is the maximum momentum also for the secondary CR 
positrons produced in the SNR. As for (2), it can be proved in the following way. First we 
can approximate the source function for secondary positrons Q 2jtt (p) as 

dlnp' 

'Pi„j l-exp(-£ e +Uil/D) 

where we assume that the spectrum of secondary positrons has the same form as that of 
primary particles except for an energy-shift factor of £ e +. Here TZ is a dimensionless constant 
which should be determined by the microphysics. Then, from Eq. ( 126)) . the escape flux of 
secondary CR positrons </> 2 (p) can be written as 

<hip) = Z(t) • $(A)PM0 + C(t), (35) 



Q2,uip) - Tlcn gas a pp K(t) exp 



"7 



(34) 
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where Z(t) and C(t) are functions of time t and do not depend on p, and $(?/) is a function 
which does not depend on time. In other words, we can rewrite 2 (p) so that the momentum 
p appears only in the combination of D p/uil. This fact shows that the peak momentum of 
2 (p) (i.e. the momentum at which d4> 2 (p)/dp = 0) is given by p 2 ,m = VmUil/D , where y m 
is the value of y at which $(?/) attains its peak value. Comparing this with Eq.flH]), we can 
say that the maximum momentum of secondary positrons decreases in the same way as that 
of primary protons, p 2 ,m oc t~ a - 



Now we can numerically calculate the escape flux of secondary positrons and its time 
evolution using the equations presented above. In Figure 3 we show the time integrated 
spectra of escaping protons and positrons. We can see that the spectral index of escaping 
secondary positrons is softer than that of escaping protons. Considering that the spectral 
index of primary electrons would be the same as primary protons, this fact shows that 
the process discussed here cannot make the rise in the positron fraction (n e +/(n e - + n e +)) 
with energy. This result is on the contrary to those presented in Blasi (2009) and Ahlers 
et al. (2009), in which they argued that the positron production by hadronic interactions 
during the acceleration of primary protons inside the SNR can reproduce the positron excess 
observed by PAMELA. 

This result can be interpreted in the following way. As shown in Fig. 1, the peak of the 
escaping positron flux is dominated by the first term of Eq. (I26[) . Using the approximation 
Eq. (J3ID, we can see that f2,o(p) oc p _7+1 for p <C Uil/D , while /2,oG°) ^ exp(— •yD p/uil) 
for p ^> u\1/Dq. Especially, at the peak momentum, the normalization factor of 2 would 
evolve with time approximately as oc ui ■ (D /u1)K(t)p2^' 1 . Considering that the spectrum 
of escaping particles should be determined by integrating Eq. ffTTj) with time, and that 
i? s h oc t 2 / 5 , ui oc t~ 3 / 5 , p 2 . m oc t~ a , and D oc t"" 1 / 5 , we can estimate the power-law index of 
the time-integrated energy spectrum of escaping positrons s esCi2 in the similar way as that 
of primary protons: 



which is obviously softer than that of escaping primary protons shown in ffl2|) . This ad- 
ditional softening, s esCj2 — s esc l = 1/a, is due to the additional factor D(p 2;rn )/u1 which 
increases with time as t . Because the momentum of secondary particles is roughly propor- 



3.2. Time Integrated Spectrum of Escaping Positrons 




(36) 
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tional to that of their parent particles, this factor is proportional to the acceleration timescale 
of protons which produce the positrons with momentum of P2, m - Especially, this timescale is 
comparable with the escape timescale, t esc ~ R^ h /D. Roughly speaking, the amount of sec- 
ondary CR particles produced in the SNR is proportional to the time during which primary 
CR particles interact with the ambient gas in the SNR. For lower energy primary protons, 
they would stay in the SNR longer (i.e. their escape timescale is longer), and so they can 
produce more secondary positrons than higher energy protons, which would stay in the SNR 
only for a shorter time and escape it into the ISM earlier. As a result, the normalization 
factor of the escaping positron flux increases faster than that of the primary proton flux, and 
therefore the time-integrated positron spectrum would be softer than the proton spectrum. 
The difference between the spectral indices of primary protons and secondary positrons is 
determined by the index a, which characterizes the evolution of the maximum momentum of 
escaping CR particles. Physically, we can say that the slower the decay of magnetic field and 
magnetohydrodynamic turbulence around the shock which determines the evolution of the 
diffusion coefficient there is, the softer the spectrum of escaping secondary particles would 
be. 



4. Discussion 

The results presented above show that the positrons which are produced during the 
acceleration of protons inside the SNR and emitted into the ISM would have a softer energy 
spectrum than those of primary protons and electrons. As is obvious from the discussions 
above, the same logic can be applied to any secondary particles which are also produced 
during the acceleration of protons in the SNRs. Hence, the emergent secondary electron 
plus positron cannot reproduce the 'excess' in the electron spectrum reported by Fermi LAT 
(Abdo et al. 2009a; Ackermann et al. 2010) because they should have a softer spectrum than 
primary electrons. In addition, the spectrum of antiprotons produced during the acceleration 
of protons and emitted from SNRs would be softer than primary protons, and so the ratio 
between them p/p would decreases with energy, which is on the contrary to the results 
presented in Blasi & Serpico (2009). 

Our results are partly consistent with those presented in previous studies such as Blasi 
(2009). In fact, in our formalism, the momentum distribution function of positrons inside 
the SNR is expressed as Eq. (1231) . and we can see that in the limit of / — > oo this is consistent 
with Eq.(6) in Blasi (2009). As this function is approximately proportional to ~ p~ J+1 when 
D(p) is linearly dependent on p, they conclude that the secondary positron (and electron) 
spectrum would be harder than the primary electron spectrum. However, this distribution 
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function is generally different from that in the ISM. For example, taking into account the 
time-dependent escape of CR particles from a SNR, even the spectrum of primary CR 
particles would be softer than that inside of the source because more and more particles 
are injected into the acceleration process as the shock sweeps the ISM (Ohira et al. 2010). 
In the case of secondary CR particles produced in a SNR, not only the increasing swept-up 
ISM gas particles injected into the shock, but also the evolution of the factor D /u^ h in 
the distribution function of secondary particles, which does not appear in that of primary 
particles, is also responsible for the softening of the energy spectrum of CR particles in the 
ISM. Considering the deceleration of the SNR shock and the decay of the magnetic field and 
turbulence around the shock, it is quite natural to assume that this factor increases with 
time, and this is just what has never been taken into account in the previous studies on this 
secondary CR production process in a SNR. As is mentioned in the previous section, because 
the lower energy positrons escape the SNR later, this factor should be larger and then the 
number flux of escaping positrons would be larger, which makes their time-integrated energy 
spectrum softer. 

We should mention the total energy of CR positrons emitted from SNRs. As the spectral 
indices of CR protons and positrons are softer than 2, most of the energy of these CR particles 
are carried by those with the lowest energy e m m, and therefore their total energy can be 
estimated as Eqr ~ £ Li n ^-^/^ £ lemm- F rom Fig- 3, where we calculate the spectra with typical 
environmental parameters (see Sec. 3), we can say that the ratio between the total energy of 
CR positrons escaping the SNR to that of secondary CR protons is < 0.06% in our case, which 
is a significant fraction of observed CR positron energy ratio to protons, ~ 0.1%. Considering 
the positron production process presented here, this ratio should be proportional to the 
product of the number density of the ISM and the escape time of the lowest energy particles, 
which is equal to the end of the Sedov phase: £?cR,e+ /Ecr, p oc nct esc (e min ) ~ nc x 10 2,5 £s- 
From Eq. fl28p . we can derive the dependence of the amount of positrons on the environmental 
parameters of a SNR as oc E s - N 1/2 n 2 / 3 ;ifthe supernova is less energetic or if the number density 
of the ISM around the supernova is larger, the amount of secondary positrons emitted from 
the SNR becomes larger. 

Our discussion may be applied also to the acceleration and escape of secondary CR nuclei 
such as boron nuclei (B) which can be produced by spallation of carbon nuclei (C) inside 
the SNR (Mertsch & Sarkar 2009). It may be worthwhile to calculate the total number 
and energy spectrum of CR borons emitted from SNRs within this scenario and give the 
theoretical lower limit to the B/C ratio in the CR observed at the Earth, which will be 
explored in the wider energy range by AMS-02. This is beyond the scope of this paper and 
left to the future work. 
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5. Summary 

In this paper we have investigated the acceleration and time-dependent escape of sec- 
ondary positrons during the acceleration of primary protons in the SNR for the first time. 
Especially we show that, if we assume that the CR particles with higher energy would escape 
earlier than lower energy particles into the ISM, the energy spectrum of secondary-produced 
particles escaping the SNR would be inevitably softer than that of primary protons. This 
implies that the rising of the positron fraction in the energy range of ~ 7 — 200GeV mea- 
sured by PAMELA and Fermi LAT cannot be accounted by the positron production during 
the acceleration of CR protons inside the SNR, and therefore it is implied that in order to 
account for this anomaly we need other primary sources of positrons with a hard spectral 
index, such as pulsars, dark matter particles, and so on. With the similar discussion, we can 
conclude that the CR electron plus positron spectrum cannot be hard enough to reproduce 
the results of Fermi LAT, and that the antiproton ratio p/p does not rise with energy within 
this scenario. These results are on the contrary to those presented in previous studies which 
have investigated this process. 

NK is grateful to K. Ioka, Y. Ohira, and T. Piran for useful discussions and their precious 
comments. This work is partially supported by an ERC advanced research grant. 
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Energy e [GeV] 

Fig. 1. — Escape flux of primary CR protons 0i(p) (thin solid line) and secondary CR 
positrons 02 (p) (thick solid line) at the time t = 30ts- The dominant components which 
make double peaks in the positron flux are also shown; 02, a (p) (dashed line) and 02,e(p) 
(dotted line). 
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Fig. 2. — Energy distribution functions of escaping CR protons (upper panel) and positrons 
(lower panel) emitted from a SNR per unit time at t — 10tg (solid lines), 30tg (dashed lines), 
100tg (dotted lines) and 300ts (dot-dashed lines). The adopted parameters are the same as 
in Fig. 1. 
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Fig. 3. — Time- integrated energy spectra of protons (black) and positrons (red) escaping a 
SNR. For convenience, the spectra are multiplied by e 2 ' 2 , and as for the positron spectrum 
the overall factor of 1500 is multiplied. The adopted parameters are the same as in Fig. 1. 



